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Abstract 

The world's marine and coastal ecosystems are under threat and single-sector management efforts have failed to address 
those threats. Scientific consensus suggests that management should evolve to focus on ecosystems and their human, 
ecological, and physical components. Estuaries are recognized globally as one of the world's most productive and most 
threatened ecosystems and many estuarine areas in British Columbia (BC) have been lost or degraded. To help prioritize 
activities and areas for regional management efforts, spatial information on human activities that adversely affect BC's 
estuaries was compiled. Using statistical analyses, estuaries were assigned to groups facing related threats that could 
benefit from similar management. The results show that estuaries in the most populated marine ecosections have the 
highest biological importance but also the highest impacts and the lowest levels of protection. This research is timely, as it 
will inform ongoing marine planning, land acquisition, and stewardship efforts in BC. 
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Introduction 

Estuaries are one of the world's most biologically productive 
ecosystems [1,2,3], providing ecosystem services valued as high as 
$22,000 US per hectare per year [4]. In British Columbia, Canada 
(BC), calculations have shown that estuarine habitats sequester a 
minimum of 180,000 tonnes of carbon annually [5] and Kelsey [6] 
estimated that 80% of all coastal wildlife use estuaries. However, 
anthropogenic threats have made estuaries one of the most 
degraded ecosystems on earth [3,7] and these threats are expected 
to increase in prevalence and severity as coastal populations 
continue to expand [7]. 

Many estuarine threats originate within the watersheds that 
drain through the estuaries [7,8]. The conversion of land for 
agriculture, forestry, or residential and industrial development 
causes habitat loss or fragmentation [7,9]. These changes in land 
use can trigger an increase in sediment, debris, nutrient, and 
pollutant levels in estuaries that smothers or entangles organisms 
or causes eutrophication, hypoxia, and anoxia [1,3,7,10,11]. 
Development that expands impervious surfaces within the 
watershed can increase water flow into estuaries, particularly 
during storm events [7], while dams impede water inputs [1,7,9]. 
As the climate warms, increased precipitation and more frequent 
storm surges are also anticipated in BC [9]. All of these 
hydrological changes can disrupt water filtration [7] and alter 
estuarine productivity and trophic structure [12,13]. 

Other threats originate in the marine environment. Develop- 
ment and industrial activities within estuaries, such as log- 
handling, aquaculture, dredging, and vessel traffic, can result in 



habitat loss or degradation, environmental contamination, shore- 
line erosion, the re-suspension of sequestered carbon, and the 
invasion of alien species [1,7,9,10,12]. Commercial and recrea- 
tional fishing can alter the composition of estuarine species 
through the preferential removal of target species [7,12]. Shoreline 
armouring to stabilize and protect coastal developments limits the 
ability of estuaries to retreat inland as the climate warms and sea 
levels rise and may therefore cause habitat and shorebird losses 
[1,9,13,14]. Higher sea temperatures may increase the severity of 
eutrophication in estuarine salt marshes [15]. Lasdy, the 
acidification of estuarine waters may inhibit the growth of shelled 
organisms [16] and increase the bioavailability of trace metals 
from pollution [17]. 

The continuing degradation of estuarine habitats has led to a 
variety of research efforts. In BC, the Pacific Estuary Conservation 
Program (PECP), a partnership of government and non-govern- 
ment organizations, has mapped the boundaries of 442 estuaries 
and designed regional scale projects to prioritize estuaries for 
conservation based on their biological importance [2]. Further 
research has found that 38% of BC estuaries had economic 
tenures within the intertidal area [18]. Merrifield et al. [1] recendy 
completed a regional assessment of threats to estuaries along the 
west coast of the United States which face similar threats to BC 
estuaries [19]. The researchers found that only 16% of surveyed 
estuaries faced few or no anthropogenic threats while 25% of 
estuaries had moderate levels of all assessed threats [1]. However, 
a systematic investigation of human activities affecting BC 
estuaries has not yet been completed. 
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In BC, the management of estuarine areas and activities 
requires the coordination of several government agencies. Subtidal 
and intertidal areas are typically public land while upland areas 
can be publicly or privately owned. The BC government regulates 
the use of public lands by issuing tenures but the ongoing activities 
of some industries with tenures, such as aquaculture, are managed 
federally. Federal jurisdiction also pertains to all organisms and 
most activities within the water column, areas designated for port 
facilities, and activities on or beneath the seabed in offshore areas 
[20]. In public areas, the BC government and several federal 
agencies have established protected areas that, to varying degrees, 
shield estuarine habitat from human activities [21]. In private 
areas, the PECP works to secure and rehabilitate important 
estuarine habitats through direct land purchases [2,22] . 

Recent reports attribute estuary degradation to a failure of 
management efforts that have focused predominandy on individ- 
ual industrial sectors [23,24,25] without examining the cumulative 
impacts of activities [24,26,27,28] or the ecosystem services on 
which the activities rely [29]. There is a growing consensus that a 
broader management focus, termed ecosystem-based management 
(EBM), is necessary to tackle the problems facing marine and 
coastal environments [23,24,26,27,30,31,32]. EBM is place-based 
and considers both the spatial variations of, and the interactions 
between, the ecological, physical, and human elements of 
ecosystems [24,31]. In Canada, Fisheries and Oceans Canada 
(DFO) has developed a framework for estuarine, coastal, and 
marine management following the principles of EBM [23]. 

To implement EBM effectively, ecosystem elements must be 
mapped and the impacts of human activities assessed 
[23,24,28,33]. This research uses geographic information systems 
(CIS) to investigate the spatial distribution of human activities that 
affect BC estuaries, incorporating threats relevant to BC that have 
not been considered previously. Cluster analysis is used to 
categorize estuaries based on the presence and magnitude of 
threat variables. Unlike in previous research [1], the spatial range 
of influence of each threat is incorporated to generate a more 
realistic spatial representation of each threat [33], the cluster 
analysis uses a non-hierarchical approach to allow the use of non- 
normal variables, and results are compared to assessments of each 
estuary's biological importance and conservation status. By 
identifying estuaries facing similar threats, this work gives context 
to previous prioritization efforts that focused solely on ecological 
factors [2], informs regional planning efforts, and helps prioritize 
limited conservation and research resources to ensure vital 
ecosystem services are sustained. 

Methodology 

Study Area 

BC's estuaries, comprising less than three percent of the 
coastline, form the study area for this research (Figure 1). I used 
estuaries and associated watersheds as the basic spatial units for 
the analysis, referred to as "estuary-watershed systems", because 
they are considered the complete geographic boundaries for 
estuarine management [8,34]. 

Estuary and Watershed Data Preparation 

In a desktop CIS, I overlaid spatial data for watersheds and 
estuaries (Table SI) to identify the watershed associated with each 
estuary (Figure 1). I confirmed that overlapping estuaries and 
watersheds shared the same primary waterway using attribute 
information found in both datasets. Estuaries that could not be 
linked to a unique watershed were removed. I added attribute 
information on importance for waterbirds [2] and carbon 



sequestration [5] to each estuary based on estuary name. I also 
overlaid the estuary-watershed systems with spatial data on 
legislated protected areas and privately conserved lands (Table 

51) to determine the count and total area of conservation 
designations within each system. Lastly, I overlaid the systems 
with BC's marine ecosections, which are planning areas with 
similar physical and biological characteristics, to determine the 
average estuarine importance and conservation values by ecosec- 
tion (Table 1). 

Threat Data Preparation 

To identify the primary anthropogenic threats facing BC 
estuaries, I reviewed the relevant literature, focusing on research 
on estuaries and estuarine habitats, such as salt marshes, seagrass 
beds, and sand flats, from BC and the Pacific coast of North 
America (Table 2). Spatial data for the primary threats (Table S2) 
were chosen based on their relevance to estuaries, spatial coverage, 
quality, resolution, and currency and used to create a compre- 
hensive spatial layer for each threat along the BC coast. While 
gaps did exist (Table S3), data were available for the majority of 
the identified threats. 

Several spatial analysis methodologies were used to compile the 
data (Table S2) for further statistical analysis. I classified land use 
data into agricultural or urban areas, combining the urban areas 
with railroads and road networks. I identified armoured shorelines 
by combining dikes and manmade shoreline segments. I selected 
the locations of freshwater diversions from a point dataset of 
obstructions. To identify areas that had been clearcut over the past 
decade, I extracted polygons with a stand age of 10 or less from a 
provincial dataset of harvest depletion. I located sources of point 
and nonpoint watershed pollution from data on pulp mills, mines, 
and areas closed due to contamination. Lastly, I calculated the 
percent change in modeled future precipitation using data on 
current (1961-1990 normal) and future (2080's CGSM A2x 
scenario) precipitation from ClimateBC [36] and the methodology 
employed by Merrifield et al. [1]. 

The compilation of marine threat data (Table S2) was guided by 
the best available data identified by the BC Marine Conservation 
Analysis (BCMCA) [37]. I obtained provincial tenure data to 
identify public areas leased to marine-related industrial activities 
and confirmed site usage with site-specific finfish aquaculture 
production data from DFO or by visually confirming the presence 
of infrastructure, such as aquaculture pens or log storage sites, 
using Google Earth 6.1 [38] satellite imagery. I merged the tenures 
to create features for shoreline development and aquaculture. I 
transformed the tenures for marinas, commercial wharves, and 
ferry terminals into points and combined them with BCMCA 
datasets representing marinas, anchorages, and tow boat reserves. 
I obtained the point locations of dredging operations and extracted 
data for large vessels from a gridded dataset of vessel traffic (Table 

52) . Smaller vessels were not well represented in this dataset so I 
used linear data on recreational boating routes to represent the 
distribution of small vessels. I combined fisheries datasets to create 
features for commercial and recreational fisheries. To retain the 
catch values available in the commercial fishing datasets, the 
values for all fisheries were converted to kg and summed. No high 
resolution datasets were available for projected ocean acidification, 
sea surface temperature, or sea level rise within BC's coastal zone 
(Table S3). 

Threats within Estuaries 

The maximum spatial range of influence of each activity was 
determined from zone values specified by Ban et al. [33] and 
Quigley and Hall [39] (Table S2). Incorporating zones more 
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Figure 1. Study area. Description: This research investigates the presence of threats within BC estuaries, which are found along the BC coastline 
from the southern border with Washington to the northern border with Alaska. The inset map shows a typical estuary and watershed. 
doi:1 0.1 371 /journal. pone.0099578.g001 



realistically represents the full spatial extent of each threat [33] and 
ensures that all relevant threats are considered, though the spatial 
footprint of associated infrastructure may not overlap directly with 
an estuary or watershed. I buffered each threat by the zone value 
to create polygons representing the larger zone of influence. 
Spatial ranges were not calculated for features with low data 
resolution [33], features where the zone of impact was not known, 
or features that already incorporated the zone of impact (Table 
S2). The spatial range was calculated for large vessel traffic, despite 
the low data resolution, due to the wide influence of this feature 



and because the original data had been fit to a coastline that did 
not incorporate the mapped estuaries in all areas. 

To determine the impact of the threats within each estuary- 
watershed system, I overlaid the buffered threats with the estuaries 
and watersheds. To standardize for variation in the size of the 
systems, I calculated the proportion of each system affected by the 
activity. For point features, I determined the density of points 
within each estuary-watershed system. For a few gridded features, 
information on level of use was available and used to better 
represent their impact spatially. Because the impact of vessel 
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Table 1. Average conservation importance, protection level, and threat status of BC estuaries, grouped by marine ecosection. 



Count of Proportion of Average 
Importance for Count of Protected Area Watershed Average Estuarine 

Marine Ecosection Waterbirds Threats Areas Protected (%) Area (ha) Area (ha) 



Continental Slope 


3.4 


6.2 


1.6 


57.9% 


4,328.3 


17.0 


Dixon Entrance 


3.2 


6.1 


1.2 


42.4% 


6,646.5 


140.0 


Hecate Strait 


3.5 


7.0 


1.3 


22.4% 


3,859.1 


55.7 


Johnstone Strait 


3.6 


8.9 


0.8 


12.8% 


5,724.1 


64.0 


Juan de Fuca Strait 


3.7 


9.1 


2.3 


11.5% 


5,019.6 


25.1 


North Coast Fjords 


3.6 


6.6 


1.3 


27.8% 


6,514.0 


192.7 


Queen Charlotte Sound 


3.5 


5.2 


1.3 


20.3% 


4,101.3 


29.4 


Queen Charlotte Strait 


3.3 


9.3 


0.9 


3.1% 


5,542.8 


125.4 


Strait of Georgia 


3.2 


11.0 


2.4 


4.2% 


7,132.2 


783.3 


Vancouver Island Shelf 


3.4 


8.9 


1.0 


13.3% 


5,255.5 


45.9 
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traffic, such as ship wake, can extend into an estuary from outside, 
I spatially joined the vessel traffic data with the estuaries and 
summed the number of vessel hours for all overlapping grid cells to 
calculate the total impact. For commercial fishing, I intersected the 
grid of fisheries catch data (Table S2) with the estuary polygons 
and summed the total catch of the overlapping grid cells 
proportional to the area of overlap to estimate the catch taken 
within each estuary. I calculated the average precipitation change 
for each estuary-watershed system using zonal statistics, which 
examined and combined the values of cells in the gridded 
precipitation dataset that fell within the estuary-watershed systems. 

I counted the number of threats affecting each estuary and 
determined the proportion of the estuary and watershed area 
covered by those activities. The average values were compared 
among estuaries in different marine ecosections (Table 1) and 
estuaries with different conservation status and biological impor- 
tance. 

Statistical Analyses 

I ran cluster analyses using Systat 1 3 [40] to group the estuaries 
based on the presence and magnitude of the threats (Table 3). 
Pairwise correlations between threats were examined to ensure the 
included threats were independent and considered equal in the 
analysis. Most of the threat variables had negatively skewed 
distributions and were log-transformed to minimize the impact of 
a few scattered higher values [28]. To determine the optimal 
number of clusters, I performed preliminary analyses with 
hierarchical clustering and the Ward method. 1 used non- 
hierarchical K-medians for the final cluster analysis because the 
method is resistant to non-normal variables and oudiers [41]. 
Euclidean distance was used to determine the distance of each site 
from the centers of each cluster. 

I mapped the results and generated spatial statistics for the 
clusters, including the average number of threats and the 
proportion of each estuary affected by each threat (Tables 3 and 
4). I also analyzed the distribution of the clusters among 
ecosections, sites with a conservation designation, and estuaries 
considered biologically important (Table 4). 



Results 

Importance and Conservation Status of Estuaries 

The average importance of estuaries for waterbirds was similar 
among all ecosections but priority carbon sequestration sites were 
only found within the Strait of Georgia and North Coast Fjords 
ecosections. Federally and provincially legislated protected areas 
overlapped 58.5% of estuary-watershed systems, predominandy in 
the North Coast Fjords ecosection. On average, 21% of estuarine 
and watershed area was covered by a legislated protected area, 
with the highest proportions in north coast ecosections. A further 
27 estuary-watershed systems, primarily on the south coast, 
overlapped a privately conserved area, with an average of 1.4% 
of the estuarine and watershed area protected. In total, 6 1 % of 
estuary-watershed systems had at least one private or legislated 
protected area. 

Threats in Estuaries 

Of 442 possible estuaries, 376 were joined to adjacent 
watersheds and included in the analysis. The literature review 
uncovered 20 threat categories and, based on the available data, 
16 threat variables were created (Table S2). A spatial dataset 
showing the threat distribution within each estuary is available 
from the author. Change in future precipitation was the only 
threat to overlap with all estuary-watershed systems but all systems 
overlapped at least two threats (Figure 2). No system overlapped all 
of the threats, though two systems in the Strait of Georgia 
ecosection were affected by 15 threats. Overall, there was an 
average of 7.9 threats per estuary-watershed system but systems 
within the Strait of Georgia ecosection averaged 1 1.0 threats and 
four other ecosections, also on the south coast, had a higher than 
average number of threats per estuary-watershed system (Table 1). 

Estuary-watershed systems with the highest importance to 
waterbirds had the highest counts of threats. Systems with a 
protected area, either legislated or privately conserved, faced an 
average of 7.6 threats. Estuary-watershed systems with less area 
protected faced a greater average number of threats and those 
facing the fewest threats were the most protected. Systems with a 
legislated protected area faced fewer average threats than those 
with a privately conserved area. 
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Figure 2. Count of threats by estuary-watershed system. Description: Estuary-watershed systems are shown classified based on the number of 
the 1 6 possible threats found within their bounds. The classification uses Jenks natural breaks which maximizes the variance between the classes. The 
number of threats faced by an estuary is greatest in the more heavily populated south coast ecosections. 
doi:1 0.1 371 /journal.pone.0099578.g002 



Cluster Analysis 

The cluster analysis revealed three optimal clusters (Figure 3). 
Cluster 3, representing estuaries with the greatest importance for 
waterbirds and carbon sequestration (Table 4), was the cluster with 
the greatest impact from most of the threat variables (Table 3). 
These 69 estuaries were found primarily in the south coast 
ecosections, dominated by the Strait of Georgia ecosection. 
Cluster 3 had the greatest average number of protected areas 
within each estuary, but the smallest average proportion of area 
protected. Cluster 2 represented 127 estuaries moderately affected 
by most threat variables and the highest impact from forestry and 
recreational fishing. The majority of estuaries in this cluster were 
found on the south coast in the less populated Vancouver Island 
Shelf ecosection. Cluster 1 faced minimal impacts from most 
threats but had a higher than average future change in 
precipitation. Cluster 1 also had the most estuaries, located 
predominantly within the North Coast Fjords ecosection. Based on 



average area conserved, these estuaries were also the most 
protected. 

Discussion 

A variety of local approaches have been important for estuary 
conservation, including the creation of legislated protected areas, 
the purchase of sites for private conservation, and the establish- 
ment of estuary management plans. A conservation designation 
was found within 61% of the estuary-watershed systems but 
systems in the north coast ecosections had greater protection than 
their southern counterparts (Table 1). The biased placement of 
legislated conservation designations in less threatened areas further 
from human settiements, and therefore less costly to protect, 
matches global trends for terrestrial protected areas [42]. Estuary- 
watershed systems with privately conserved areas faced a greater 
number of threats than those with legislated protected areas. As 
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Figure 3. Distribution of K-medians cluster results. Description: Estuary-watershed systems are shown distributed into clusters that resulted 
from the K-medians cluster analysis based on the presence of 1 6 threat variables. Similar to the count of threats, the clusters highlight a north-south 
gradient, with highly threatened estuaries primarily observed in the south coast ecosections and minimally threatened estuaries predominantly 
observed in the north coast ecosections. 
doi:1 0.1 371 /journal.pone.0099578.g003 



Brooks et al. [42] note, people most value "what is close to them". 
Privately protected areas, often negotiated by the PECP with 
landowners, are found primarily along the south coast, where the 
greatest concentration of humans, and industrial activities, is found 
[10,17]. 

In 2002, a review of estuary management plans in BC identified 
only nine completed plans, predominantly for estuaries on the 
south coast [43] . Since that time, planning has been completed for 
two additional estuaries but some of these plans, particularly those 
with limited stakeholder support, have taken up to 20 years to 
finalize [43] and some have not been fully implemented (K. 
Moore, Environment Canada, personal communication). With 
uneven estuarine protection along the coast and management 
plans for less than 2% of mapped estuaries, it is apparent that local 
efforts alone will not be sufficient to conserve estuarine habitats 
effectively into the future. 



To ensure the maintenance of biodiversity, a regional approach 
based on the principles of EBM offers a more comprehensive, 
efficient, and cost-effective alternative to site-specific protection 
and planning [30,44]. The PECP's efforts to map and rank 
estuaries began the work of looking at BC's estuaries through a 
regional lens. However, the PECP lacked the comprehensive 
anthropogenic data necessary to complete the ranking [2]. 
Identifying clusters of estuaries that face similar threats can 
highlight priority activities and areas where coordinated manage- 
ment could achieve a higher return on investment. 

Highly Threatened Estuaries 

Areas with high diversity or importance and high cumulative 
impacts are considered conservation priorities to protect species 
and habitats that would otherwise be lost to habitat modification 
[28,42]. Identified as highly important for waterbirds and carbon 
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sequestration, the highly threatened estuaries in Cluster 3 qualify 
as conservation priorities (Table 4). These estuaries were located 
predominantly within the heavily populated south coast ecosec- 
tions (Figure 3), which correlates with results from previous 
research [18,33] as well as the simpler count of threats by estuary 
(Table 1 and Figure 2). 

Protected areas that prohibit human activities are a prominent 
means of conserving estuaries [23] but the systems in Cluster 3 had 
the lowest average proportion of area protected (Table 4). Many of 
those estuaries are located near major urban centers, where the 
financial costs and logistical challenges associated with establishing 
and ensuring compliance of conservation designations are higher 
[45]. Protecting estuaries in highly populated areas will require 
stakeholder engagement, a component of EBM, to ensure that all 
relevant ecosystems and activities are considered, that priorities 
are established collectively, and that efforts are not duplicated 
[30,46]. However, only the conservation of representative sites 
across the coast will ensure the continued viability of estuarine 
ecosystems and help Canada meet its international commitment to 
preserve 10% of each ecological region in the marine and coastal 
realm [46]. 

Moderately Threatened Estuaries 

Cluster 2 faced a moderate level of most threats (Table 3) and 
contained estuaries located primarily in the less populated 
ecosections of the south coast (Figure 3). These systems had the 
greatest impact from recreational fishing and forestry, suggesting 
that coordinating estuarine management with these two sectors 
might alleviate significant stresses. Regional planning is a means of 
engaging industries whose activities adversely affect an ecosystem 
from afar [47], such as the terrestrial activities like agriculture, 
forestry, and watershed development that affect estuaries in 
Clusters 2 and 3. Smith el al. [47] suggest that a system of local 
plans for highly used areas nested within regional plans might offer 
the most efficient solution. This suggestion meshes with recom- 
mendations from the Royal Society of Canada [48], who proposed 
that planning for DFO's large ocean management areas should be 
tied to planning for the smaller coastal management areas because 
the areas are physically and biologically integrated. 

Minimally Threatened Estuaries 

While some prioritize action in highly threatened areas, there 
are also calls to minimize conflict by proactively protecting areas 
required by few industries [49,50]. Over 40% of the estuary- 
watershed systems assessed, represented by Cluster 1, were 
minimally affected (Table 3). This result is encouraging, suggesting 
that many of BC's estuaries are relatively pristine. In comparison, 
research done for estuaries along the more heavily populated west 
coast of the United States found that only 16% were minimally 
affected [1], The systems of Cluster 1 were located primarily in the 
lightly populated north coast ecosections (Figure 3). These sites 
faced a slighdy greater average increase in future precipitation but 
because of the relatively low levels of other anthropogenic threats, 
the estuaries may be best poised to respond to the corresponding 
increase in water flow [51]. 

While the estuaries of Cluster 1 were the least affected by the 
combination of all threats, they did show moderate impacts from a 
few threats (Table 3). These estuaries are primarily located on the 
north coast, away from human setdements. However, small 
human populations do not guarantee minimal anthropogenic 
impacts [28] and some threats, such as forestry, shipping, and 
climate change, can influence even remote estuaries. Threats that 
impact all of the clusters are examples of regional issues that are 
likely managed most efficiently at a broad scale [1 1]. 



Limitations in Data and Analysis 

Cluster analysis is an inherently subjective statistical methodol- 
ogy relying on researcher judgment to gauge the integrity of the 
results. A strong theoretical foundation for the analysis helps 
exclude irrelevant variables that might distort results [41] and the 
extensive literature review of anthropogenic threats (Table 2) 
ensured only well-researched variables were incorporated. How- 
ever, the lack of fine scale or historical spatial data on marine 
activities and climate change is a noted impediment to conserva- 
tion and research efforts [17,44], 

Two specific data gaps should be noted. Observations of 
invasive species, one of the most important stresses in estuarine 
ecosystems [9,19], have not been catalogued across the BC coast 
(D. Buffett, Ducks Unlimited Canada, personal communication). 
Instead, two vectors of invasive species, aquaculture and large 
vessel traffic [10], were used and the results can help identify 
highly disturbed estuaries most susceptible to invasion [8] . Perhaps 
the most important data gap is projected sea level rise. The 
shoreline armouring feature represents estuaries with barriers 
preventing adaptation to higher sea levels but does not consider 
spatial variations in elevation. The vulnerability of BC shorelines 
to sea level rise has been approximated by the BC Government 
using coarse elevation data and the physical parameters of 
shoreline segments [17]. However, better projections of the 
impact of SLR on estuaries require high resolution elevation data 
that are not currendy available for much of the coast [17]. 

Limitations in Implementation 

Following the approach taken by Merrifield et al. [1], this 
analysis incorporated threats found throughout estuary-watershed 
systems. Estuaries link terrestrial watersheds with marine systems 
and the influences of the two realms both enhance [52] and 
threaten [3,7] the ecosystem services provided by estuaries. 
However, while terrestrial, coastal, and marine environments are 
integrated, they are often not governed harmoniously. Considering 
the entire land to sea continuum is a more comprehensive 
approach to conservation and management advocated by both 
government agencies [23,27] and researchers [7,8,34,47,53], 
despite the challenges associated with the large area and number 
of stakeholders involved. A failure to utilize EBM and explicidy 
consider threats originating across the land-sea continuum hinders 
effective estuary management and can undercut the success of 
protected areas [30]. Conservation efforts require the coordinated 
involvement of all relevant government bodies to facilitate 
successful integrated and adaptive planning [7,21,23,34]. This 
will be particularly relevant for the estuaries in Clusters 2 and 3, 
which face a suite of significant marine and terrestrial threats. In 
BC, improved coordination may be underway as the BC 
government and First Nations are developing an EBM plan that 
considers the cumulative impacts of human activities (www. 
mappocean.org) and federal agencies are collaborating on a 
network of marine protected areas [46] . 

Priorities for Management and Conservation 

Land purchases such as those coordinated by the PECP play an 
integral role in the conservation of estuaries. To a certain extent, 
these purchases will likely always be done in an ad hoc manner, 
taking advantage of the generosity of landowners adjacent to 
estuaries. However, planning at a regional scale using biological 
and socioeconomic data can help ensure that management and 
conservation efforts are proactive, scientifically defensible, maxi- 
mize available resources, and minimize impacts on local resource 
users [2,25,54]. A comprehensive network approach to protected 
areas can also help guarantee that protected areas fully realize 
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their benefits [50]. To conserve estuarine ecosystem services, a 
network of protected areas must ensure that contributing sites 
include replications across the coast, maintain connectivity, and 
are adequately sized and protected [46]. 

Estuarine conservation and management would also benefit 
from further research into the cumulative impacts of estuarine 
threats. This analysis determined the specific threats facing BC's 
estuaries but more information is required to understand how 
threats interact, in particular for the highly threatened estuaries in 
Clusters 2 and 3. Ban et al. [33] began this work for BG by 
extrapolating expert assessments from California on the vulnera- 
bility of marine habitats to different activities [35]. However, 
information on the level of impact of each threat is also required 
and is currently only available for a few of the regional scale 
datasets, a shortcoming identified in other jurisdictions [1]. The 
BCMCA has worked with industry representatives in the marine 
environment to validate the accuracy of spatial data [37] and 
continuing that work to generate a better representation of site 
usage and impact is an important next step. 

Austin et al. [9] declared climate change to be the "foremost 
threat to biodiversity" in BC but a lack of spatial data (Table S3) 
make it challenging to incorporate the relevant threats. The BC 
estuaries identified as priorities for carbon sequestration [5] were 
shown to be highly threatened (Table 4). Protecting pristine 
carbon sequestration sites and restoring degraded areas through 
an EBM process has been referred to as "one of the strongest win- 
win mitigation efforts known today" [32] because safeguarding 
carbon sequestration requires maintaining estuarine structure and 
function, facilitating the continued provision of the other 
ecosystem services that estuaries provide. Emerging information 
on the capacity of estuaries to sequester carbon, along with better 
accounting methodologies, may bring renewed impetus to improve 
and fund coastal management from both the public and private 
sectors [22,32,42,55]. 

Conclusions 

Estuaries are an important component of BC's biodiversity and 
the ecosystem services they provide are well established. However, 
estuary degradation continues, in part due to a historical lack of 
integrated management and the investment required to develop 
and implement individual estuary management plans. Regional 
planning based on EBM that involves stakeholders across the 
entire land-sea continuum is required to prioritize sites for 
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